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ABSTRACT 


Three-dimensional fully developed turbulent mixed convection heat transfer of Al,03/synthetic oil 
nanofluid in a trough collector tube with a non-uniform heat flux was numerically studied. The effect of 
AL203 particle concentration in the synthetic oil on the rate of heat transfer from the absorber tube was 
also investigated. Various nanoparticle concentrations ( < 5% in volume) at the operational temperatures 
of 300 K, 400 K, and 500 K were used in the current study. The heat flux distribution on the outer surface 
of the absorber tube was non-uniform in the circumferential direction but uniform in the axial direction. 
The heat flux in the circumferential direction was obtained using the Monte Carlo ray tracing technique. 
Three-dimensional Navier-Stokes mass, momentum and energy equations were solved using the FLUENT 
software. A second order discretization was used for the convective and diffusive terms and the SIMPLE 
scheme was applied for the velocity-pressure coupling. The outlet temperature of the absorber tube 


Benen from the current numerical simulations was compared with existing experimental data. The absolute 
errors were found to be less than 3.8 °C. The current numerical results show that the convection heat 
transfer coefficient has a direct dependency on the volumetric concentration of nanoparticles in the fluid. 
In addition, the heat transfer enhancement due to the nanoparticles in the fluid reduces as the absorber 
operational temperature is increased. 

© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


90% of the world's energy demands. The fossil fuels are non- 
renewable and limited in supplying energy indefinitely. These 


The global energy demand is increasing in an accelerating rate. energy resources will eventually be depleted. In addition, combus- 
Fossil fuels (oil, coal and natural gas) are currently supplying over tion of fossil fuels makes carbon dioxide emission, which is one of 


the greenhouse gases contributing to global warming. Use of 
renewable energy resources, on the other hand, can lead to reduction 
in fossil fuel consumption and in turn protection of the environment. 


* Corresponding author. Tel.: +98 9121947510; fax: +98 21 88617087. 


E-mail address: akasa@ut.ac.ir (A.B. Kasaeian). Solar energy is one of the largest renewable energy resources. 


http://dx.doi.org/10.1016/j.rser.2014.02.028 
1364-0321 © 2014 Elsevier Ltd. All rights reserved. 


T. Sokhansefat et al. / Renewable and Sustainable Energy Reviews 33 (2014) 636-644 637 


Nomenclature 


(ar Specific heat capacity (J/kg K) 

g Gravitational acceleration 

h Heat transfer coefficient (W/m? K) 
K incident angle modifier 


K Thermal conductivity (W/m K) 
L Length of tube (m) 

p Pressure (Pa) 

q’ Heat flux (W/m?) 

Dg Diameter of absorber tube (m) 
Re Reynolds number 

T Temperature (K) 


The greatest advantage of the solar energy over alternative forms of 
fossil energy resources is that it is clean and abundant. 

Although the costs of solar energy have gone down and 
continue to fall, the levelized costs of solar energy are still much 
higher than conventional energy. The levelized cost of solar 
concentrated power (CSP) (~ 194-336 USD/MW h) is four times 
than that of supercritical coal without carbon capture and storage 
(~43-88 USD/MW h) [1]. Cost is one of the major factors inhibit- 
ing development of trough collectors. Increasing the efficiency of 
the solar collector systems can in part alleviate this problem. Many 
attempts have been made to improve the performance of these 
systems [2-5]. 

Heat transfer can be enhanced by increasing the thermal 
conductivity of the heat transfer fluid (HTF). Thermal conductivity 
of metallic particles, metallic oxides and nanotubes is relatively 
higher than that of liquids. Addition of fine particles into heat 
transfer fluids (thus forming nanofluids) can significantly increase 
the heat transfer rate [6-8]. It is obvious that the researches on 
the applications of nanofluids have been popularized during the 
recent years. Many studies have been carried out on the effect of 
nanofluids and nanofluid size on heat transfer [9-10]. A number of 
investigations show the improvement of convection heat transfer 
using nanofluid. Choi et al. [11] studied engine oil/carbon nano- 
tube nanofluids. They observed 160% enhancement of thermal 
conductivity using 1% volume carbon nanotube. The concentration 
of nanoparticles has high effects on the fluid convection heat 
transfer coefficient [12]. Akbari and Behzadmehr [13] investigated 
the mixed convection heat transfer of a water/Al,03 nanofluid 
system in a horizontal tube with uniform heat flux. Their results 
demonstrate that the convective heat transfer coefficient is 
increased from 9% to 15% using 2% and 4% volume Al0O3 nano- 
particles, respectively. Duangthongsuk and Wongwises [14] stu- 
died a system of water nanofluid consisting 0.2% volume TiO, 
nanoparticles which caused 6-11% enhancement of heat transfer 
coefficient. 

Teng et al. [15] surveyed the changes in heat transfer of Al203- 
water nanofluid at different diameter sizes of nanoparticles. They 
declared better thermal conductivity in the smaller nanoparticle 
diameters. In another study, Heris et al. [16] compared the heat 
transfer performance of Al,O3/water with CuO/water nanofluids. 
They observed that the obtained heat transfer enhancement with 
Al,03 nanoparticles is higher than that of CuO nanoparticles. As an 
example, they stated that at a Peclet number of 5000, the heat 
transfer coefficient enhancement ratio was 1.29 for 2.5% volume of 
Al,03 nanoparticles whereas it was 1.23 for CuO/water nanofluid 
for the same particle volume fraction. Bianco et al. | 17] considered 
a laminar flow of Al,O3/water nanofluid. Under constant wall heat 
flux boundary conditions, they analyzed the heat transfer of 


AT Average fluid temperature above ambient air tempera- 
ture(K) 

u,U,w Dimensionless velocity (m/s) 

0 Angular coordinate 

r Radial direction 

H Dynamic viscosity (kg/m s) 

h Particle volume fraction 

p Density (kg/m?) 

a Absorber 

bf Base fluid 

f Focal length (m) 

eff Effective 

S Solid particle 

W Wall 


nanofluid using both single phase and two-phase approaches. 
They noted that the difference between the results of single phase 
and two-phase approaches is small, especially when the tempera- 
ture dependency of the thermo-physical properties is taken into 
account. 

Using nanofluid as a working fluid in solar systems is a novel 
approach to increase the efficiency of solar systems. Lenert and 
Wang [18] studied the influence of different variations in nano- 
fluid volumetric receivers both theoretically and experimentally. 
In their theoretical part, a one-dimensional transient heat transfer 
model was supposed and the enhancement of receiver efficiency 
with augmentation of nanofluid height (H) and incident solar 
flux was proved. In the experimental part, carbon coated cobalt 
nanoparticles were added to Therminol™ VP-1 in a liquid-based 
volumetric receiver. They illustrated that the efficiency of nano- 
fluid volumetric receivers is increased with increasing solar 
concentration and nanofluid height. Also, they stated that the 
receiver-side efficiencies are predicted to exceed 35% when nano- 
fluid volumetric receivers are coupled to a power cycle and 
optimized with respect to the optical thickness and solar exposure 
time. Lu et al. [19] evaluated the influence of water-based CuO 
nanofluids upon an open thermo syphon utilizing in a high- 
temperature evacuated tabular solar collector (HTC). Compared 
to water, nanofluid improved the thermal performance of the 
evaporator; also a 30% enhancement was observed for the evapor- 
ating heat transfer coefficient. 

Yousefi et al. [20] investigated the effect of Al203-H20 nano- 
fluid on the efficiency of a flat-plate solar collector experimentally. 
In comparison with pure water, they declared 28.3% enhancement 
of efficiency for 0.2 wt% nanoparticle concentration. Also they 
showed that utilizing MWCNT-H20 nanofluid enhanced the effi- 
ciency of a flat-plate solar collector [21]. Khullar et al. [22] studied 
the nanofluid usage in a concentrating parabolic solar collector 
theoretically and compared the results with the experimental data 
of the conventional concentrating parabolic solar collectors. Their 
results demonstrate 5-10% higher efficiency as compared to the 
conventional models. Kasaeian et al. [23] investigated the heat 
transfer enhancement for Al,03/synthetic oil nanofluid in a para- 
bolic trough collector tube numerically. Their results show that the 
heat transfer coefficient is increased as the concentration of the 
nanoparticles in the base fluid is increased. Otanicar et al. [24] 
examined the effect of different nanofluids (carbon nanotubes, 
graphite and silver) on the performance of a direct absorption 
collector. The efficiency rose with the enhancement of particles 
concentration but, after a volume fraction of 5%, the efficiency 
diminished slightly. Since the nanofluid receiver can potentially be 
either useful in a direct steam generation system (DSG), some 
works are reported in the literature [25]. 
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Table 1 
Economic comparisons for conventional and nanofluid-based solar collectors [26]. 


Parameter Conventional solar Nanofluid solar 
collector ($) collector (S) 

Capital cost 

Independent costs 200 200 

Area based costs 397.8 327.8 

Nanoparticles 188.79 

Total capital 597.8 716.59 
(one time cost) 

Total maintenance 96.23 115.35 
(for 15 year life) 

Total cost 694.03 831.94 

Electricity cost 270.13 278.95 
saving per year 

Years until electricity 2.57 2.98 
saving =cost 

Natural gas cost 80.37 83.02 
saving per year 

Years until natural 8.64 10.02 
gas saving=cost 

Electricity price 

November-March 0.08 0.08 
(per kW h) 

May-October (per kW h) 0.09 0.09 

Daily service charge 0.25 0.25 

Gas price rate (per term) 0.74 0.74 

Monthly service charge 9.70 9.70 


The use of nanofluids in solar systems is desired to be evaluated 
in the economic point of view. Otanicar and Golden [26] compared 
the economic features of nanofluid-based solar collectors with the 
conventional types. The study was based on life cycle assessment 
which is a capable methodology to evaluate the economic impacts 
of products. As it is shown in Table 1, the capital and maintenance 
costs are $120 and $20 higher for nanofluid-based collectors 
respectively for the year life time. The payback period is less for 
the conventional collectors, but with assuming 15 years life time 
according to the better performance of nanofluid-based collectors, 
the life cycle savings would be nearly the same. Taylor et al. [25] 
evaluated the utilization of nanofluid receivers for a solar thermal 
power plant theoretically. They appraised the amount of nanopar- 
ticles that would be needed for a solar thermal power plant (3 kg 
per each 1 MWe) and with assuming the price of nanoparticles 
~$1000/kg, the total capital investment would be $5/W where the 
cost rise is less than 0.1% of the total capital investments. They 
proposed two conceptual designs for nanofluid receivers that 
substitute the conventional models in the solar plant which are 
cheaper than the conventional ceramic receivers. They used 
graphite/therminol® VP-1 nanofluid in the model with a volume 
fraction less than 0.001%. The conservative calculations exhibited 
that for a 100 MWe power tower solar plant operating in Tucson- 
Arizona, more than about $3.5 million was obtained in the yearly 
revenue and the payback time of the plant reduces about two 
years. 

The objective of the current research is to model a trough 
collector absorber tube and study the effects of nanoparticle 
concentration on the mixed convection heat transfer rate of the 
nano-fluid in a fully developed turbulent flow. Most of the models, 
in the literature, for simulating the flow in a trough collector 
absorber tube assume that the solar flux is uniform and many 
correlations in the models are based on a uniform temperature. 
Because of the non-uniform nature of the solar flux on the outer 
surface of the inner absorber tube, the temperature distribution 
of the heat transfer fluid is, in fact, asymmetric. Relatively few 
investigations have been reported on simulating non-uniform heat 
flux in absorber tubes [27,28]. Al,O3/synthetic oil was used as the 
heat transfer fluid. The fluid properties are strongly dependent on 


the absorber operational temperature. The heat transfer coeffi- 
cients for various volumetric concentrations of nanoparticles in 
the heat transfer fluid are presented along the axial and circum- 
ferential directions for a number of absorber tube operational 
temperatures. 


2. Fluid properties 


Thermo-physical properties of the base fluid (synthetic oil) and 
the nanoparticles (Al203) including density, viscosity, thermal 
conductivity, and specific heat are strongly dependent on the 
operational temperature. The variation of density of the base fluid 
with temperature was extracted from Ref. [29] and plotted in Fig. 1(a). 
The distributions of viscosity, thermal conductivity, and specific heat 
of the base fluid as a function of temperature were also obtained from 
[27] and used in the current study (Fig. 1(b) and (c)). Also, the 
properties of Al,O3 were taken from Ref. [30]. Addition of nanopar- 
ticles remarkably alters the thermo-physical properties of the heat 
transfer fluid. These properties were calculated for various concentra- 
tions of nanoparticles in the fluid and plotted in the same figure for 
compassion. Table A1 of Appendix A provides the numerical values of 
these properties for the operational temperatures of 300 K, 400 K, 
and 500 K. The thermo-physical properties of the base fluid, nano- 
particles and nanofluid as a function of temperature are provided in 
Appendix A. 


3. Numerical modeling 


The flow of heat transfer fluid through an LS-2 parabolic trough 
collector module, tested on the AZTRAK rotating test platform at 
SNL by Dudley et al. [31], was numerically simulated. A schematic 
diagram of the collector and absorber tube is shown in Fig. 2 and 
the geometric characteristics of the collector are shown in Table 2. 
The absorber tube is 7.8 m (L) long and 0.07 m (D,) in diameter. 
Al,03/synthetic oil was used as the working fluid in the collector 
field and the nanofluid flow was modeled as single phase. Since 
the Re is higher than 2500 and Richardson number (Ri= Gr/Re?) is 
between 0.1 and 10, the flow is assumed to be fully developed 
turbulent mixed convection. Gravitational acceleration was also 
taken into account. The operational temperatures are 300 K, 400 K, 
and 500 K and the pressure is assumed 20 psig. The fluid motion 
governing equations were solved using the FLUENT software. The 
governing equations are provided in Appendix A. 

A k-e formulation was used to model the turbulence in 
the flow. In the k-e model, the turbulence characteristics are 
described based on the turbulent kinetic energy (TKE) and viscous 
dissipation of turbulent kinetic energy. The FLUENT software was 
used to numerically simulate the flow of the nano-fluid in the 
absorber tube. This software is based on finite volume method. 
To couple the velocity field and pressure, the SIMPLE model was 
employed. In this method, the pressure correction equation is 
solved in a number of iterations, and then the velocity is corrected 
until the continuity equation is satisfied in the computational 
domain. For the convective and diffusive terms, a second order 
discretization method was used. The flow was assumed three- 
dimensional, steady state, and single-phase. To insure that the 
results were time independent, the average heat transfer coeffi- 
cient was monitored with the size of the grid. 

The solar energy flux distribution on the outer wall of the inner 
absorber tube was obtained using the Monte Carlo Ray Trace 
(MCRT) technique developed by He et al. [27]. The heat flux 
distribution of He et al. (Fig. 3) was used in the current study. In 
the numerical simulations, the boundary condition on the absor- 
ber tube wall was set to this non-uniform heat flux distribution 


T. Sokhansefat et al. / Renewable and Sustainable Energy Reviews 33 (2014) 636-644 639 


a 
1200 T 
Base Fluid (BF) 
>= 0.01 
1100 = 0.03 
“Ee = 005 
“ob 
% 
E 1000 
n 
5 
Q 
900 
800 $ ; $ 
300 400 500 600 700 
T&) 
Cc 
0.18 
Base Fluid (BF) 
0.16 o=0.01 
= 0.03 


0.14 = 005 
0.12 
0.1 


0.08 


Thermal Conductivity (W/m.K) 


0.06 , . . 
300 350 400 450 500 550 600 


T®) 


b 


——— Base Fluid (BF) 
$=001 
$= 0.03 
= 0.05 


Viscosity (mPa.s) 


a 


2.2 
S a 
[eto] 
= 
% 
g 18 
T 
2 
3 Base Fluid (BE) 
& 1.6 = 0.01 
N 
$= 0.03 
= 0.05 
1.4 , , 
300 350 400 450 500 550 600 


TK) 


Fig. 1. (a) Density, (b) viscosity, (c) thermal conductivity, and (d) specific heat of the base fluid as well as the nanofluid for various particle concentrations at different 


operational temperatures. 


Parabola 


Fig. 2. Schematic diagram of the parabolic trough collector and absorber tube. Red 
lines represent rays of light from the Sun. (For interpretation of the references to 
color in this figure caption, the reader is referred to the web version of this article.) 


Table 2 
Specifications of the LS-2 parabolic trough collector [31]. 


System Type Parameter Value 
Parabolic trough reflector Zp%X, 7.8 m*5 m 
f 1.84 m 
Receiver Tg 0.0575 
Ta 0.035 


using a user-defined function in the FLUENT software. In addition, 
the boundary conditions at the inlet (z=0) were as 


v=0, THT 


W=Wo, u 


where u, v, and w are the velocity components in x, y, and z 
directions, respectively, and T is the temperature of the fluid. The 
flow at the outlet of the tube was assumed fully developed and an 
outflow condition was used in the simulations, the fully developed 


qi W-m? 


-90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 
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Fig. 3. Heat flux distribution of the absorber tube as a function of angle. Graph was 
reproduced from [27]. 
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Fig. 4. Axial velocity profiles at different axial positions at Re=26 000. 


entrance length is demonstrated in Appendix A. Also, Fig. 4 shows 
the axial velocity profiles at six different cross sections for a given 
Re (260 00). We can observe that the fully developed region occurs 
in Z/D=30. 
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Furthermore, the following boundary conditions were used at 
the fluid-solid interface: 


oT 
Keff Op 


u=v=w=0, q’= 


The MCRT method is dependent on the geometric (r, Z,) as well 
as the grid (Ne, Nz) parameters. The MCRT method requires a grid 
density for the absorber tube of N, > 68 and N; > 320. A structured 
non-uniform mesh was created to discretize the computational 
domain. It is finer near the wall where the velocity and tempera- 
ture gradients are larger. The standard method to test for grid 
independence is to increase the resolution and repeat the simula- 
tion. If the results do not change appreciably, the original grid is 
probably adequate and the results are not dependent on the girds 
size. If, on the other hand, there are significant differences 
between the two solutions, the original grid is likely of inadequate 
resolution. In such a case, an even finer grid should be tried until 
the grid is adequately resolved. A number of grid sizes were 
studied to insure that the accuracy of the results is independent of 
the grid size. The mean heat transfer coefficient was monitored for 
two cases: (1) T=500 K and g=0 and (2) T=500 K and m=0.05. 
The results of the grid refinement study are presented in Table 3. 
Increasing the grid numbers does not significantly change the mean 
heat transfer coefficient. The final mesh for the current study 
consisted of 68 (N<), 320 (Nz), and 20 (N,) nodes in the circumfer- 
ential, axial and radial directions, respectively. 


Table 3 

Average heat transfer coefficients at 500 K for various grid sizes. 
Grid (NéN#N-) Mean heat transfer 

coefficient at ¢=0.05 


Mean heat transfer 
coefficient at ¢=0 


68-20-320 137.99 146.53 

68-30-400 137.75 147.01 

80-20-480 138.01 146.55 
Table 4 


The experimental data from Dudley and comparisons with our simulation result. 


4. Validation 


Two available experimental cases were used to validate the 
current numerical study. 


4.1. LS-2 solar collector 


Three typical states for a bare tube from Dudley et al.'s report 
[31] were used for the purpose of validation. The results are 
presented in Table 4. Dudley et al. derived the performance 
correlation between the collector efficiency and the working fluid 
temperature as following [31]: 


2 
n = K*(73.3 — 0.007276%AT) + (496% (E) — 0.0691% Ge ) 
b b 


In the equation, K is the incident angle modifier and may be 
assumed to be one (K=1). 

The heat flux distribution on the outer surface of absorber tube 
is shown in Fig. 5(a). The temperature distribution on the outer 
surface of absorber tube is also shown in Fig. 5(b). Fig. 5(b) shows 
the wall temperature distributions of the absorber tube based on 
the test condition state 1. Fig. 6 shows the temperature distribu- 
tion in the axial direction with 0=0°, 90°, 180°, 270° and 310°. The 
wall temperature increases with the increase of z along axial 
direction, and the temperature at 310° is much larger than the 
other angular direction. The reason is that, 0=310° is on the side 
with high concentrated solar radiation, while @=90° is on the side 
without concentrated solar radiation. According to the simulated 
outlet temperature and collector efficiency based on the test 
results of Dudley (Table 4), the absolute error is less than 3.8 °C 
which has good agreement with the obtained simulation results. 
This proves that the models and methods used in the present 
study are feasible and the numerical results are reliable 


To,test (°C) To,model (°C) E (%) (test) (%) Himodel) (%) 


933.7 0.66 4501 212 102.2 


124 126.1 1.69 72.51 72.27 
173.3 176.9 2.1 70.9 71.5 


316.9 68.98 66.9 


: 319.32 0.76 


Test Ip (W/m?) tn (kg/s) Reynolds y PE la tt Tin (°C) 
State 1 
State 2 968.2 0.65 7653 22.4 151 
State 3 937.9 0.61 28560 28.8 297.8 
a 
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Fig. 5. (a) Heat flux and (b) temperature distributions on the outer surface of the absorber tube. 
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4.2. Al203 nanofluid flow through a horizontal tube 


Heat transfer coefficient and Nusselt number were obtained and 
compared with the experimental results of Esmaeilzadeh et al. [32] 
in a horizontal tube. Fig. 7(a) and (b) shows the both the results of 
the current numerical study and those of Esmaeilzadeh et al. [32]. 
The axial evolution of Nusselt number and heat transfer coefficient 
are in good agreement with the corresponding experimental results. 
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Fig. 6. Mean temperature distributions along the axial direction (z) on the outer 
surface of absorber tube at various angular locations at Re=17 361. 
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5. Results and discussion 


The study of heat transfer was carried out on a trough collector 
absorber tube with 7.8 m length and 0.07 m diameter. The tube is 
consumed to be evacuated inside the space between the absorber 
metal tube and the glass tube. The vacuum level in the annular 
space between the inner absorber tube and the outer glass cover 
tube is so high that the convective heat losses could be essentially 
eliminated, but the thermal radiation heat transfer between the 
outer wall of absorber tube and the inner surface of glass cover 
tube is taken into account. Also, the heat transfer mechanism 
inside the tube is as convection and the inlet fluid velocity is 
assumed to be 4 m/s. The convection heat transfer of nanofluid 
may be increased compare to the base fluid. The improved thermal 
conductivity of nanofluid is the most important factor for improv- 
ing heat transfer in parabolic trough collectors. Actually, nanofluid 
reduces the thermal resistance at interfaces and minimizes the 
temperature difference between the absorber and the heat trans- 
fer fluid. The effects of the particle concentration on the flow fields 
and heat transfer coefficient along the absorber tube are pre- 
sented. Convection heat transfer coefficients at various axial 
locations for a given Reynolds number and five volumetric con- 
centrations of Al20; particles are shown in Fig. 8. Because of the 
non-uniform heat flux nature, a number of graphs are shown for 
heat transfer coefficients (Fig. 8). As shown, the heat transfer 
coefficient ranges from 100 to 260 and 60 to 160 at the operational 
temperatures of 300K and 500K, respectively, indicating the 
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Fig. 7. Comparison of (a) the heat transfer coefficients, (b) the Nusselt number along the axial direction from the current numerical study with the existing experimental 
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Fig. 8. Convective heat transfer coefficient vs. axial distance for various particle concentrations and angular locations for (a) Re=3900 (T=300 K) and (b) Re=26 800 


(T=500 K). 
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Fig. 9. Mean convective heat transfer coefficient vs. particle concentration at the 
operational temperatures of 300, 400 and 500 K. 
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Fig. 10. Percentage increase in mean convection heat transfer coefficient due to the 
addition of nanoparticles to the base fluid vs. operational temperature for various 
particle volume concentrations. 


decrease of heat transfer coefficient with the increase in opera- 
tional temperature. 

The heat transfer coefficient sharply drops and then levels 
off to its fully developed value. As the nanoparticle concentration 
increases, the heat transfer coefficient rises. The coefficient is more 
affected by the particle concentration at the operation tempera- 
ture of 300 K than 500 K. However, the heat transfer coefficient is 
reduced as the operational temperature is increased from 300 K to 
500 K. The average heat transfer coefficients are lower for the 
higher temperatures (Fig. 9). As shown in this figure, the average 
convection heat transfer coefficient increases by about 14% by 
increasing the nanoparticle volumetric concentration by 5% at the 
operational temperature of 300 K. Also, the effect of nanofluid 
concentrations on average heat transfer coefficient at the opera- 
tional temperature of 400 K increases by 8.6% with the addition of 
5% (in volume) of Al2O3 nanoparticles. This effect was also studied 
at 500 K and the results show that the average convective heat 
transfer coefficient increases 6% by adding 5% (in volume) of Al203 
nanoparticles. 

From Fig. 1, we can observe that when the operational 
temperature is 300K the value of „ increases to 8.38e—4 and 
the value of k increase to 0.028 by adding 5vol% of Al203 
nanoparticles. While at the operational temperature of 400 K, 
the value of u increases to 2.7e—4 and the value of k increase to 
0.023 by adding 5 vol% of Al,O3 nanoparticles. Also at the opera- 
tional temperature of 500 K the value of „u increases to 9.6e—5 
and the value of k increases to 0.019 by adding 5 vol% of Al203 
nanoparticles. The other reason for the lower increase in heat 
transfer coefficient at higher temperatures can be associated to the 


decrease in viscosity and thermal conductivity change rates at 
higher temperatures [23]. The percentage increase in heat transfer 
at various temperatures is shown in Fig. 10. According to this 
figure, the percentage of heat transfer increase has a descending 
trend with temperature rising. 


6. Conclusions 


The thermal behavior of a nanofluid consisting of synthetic oil 
and Al,03 throughout a parabolic trough collector absorber tube 
(L=7.8 and D=0.07 m) has been simulated and a turbulent mixed 
convection heat transfer with non-uniform heat flux has been 
investigated. The results of the current study illustrate that the 
heat transfer coefficient of the working fluid in an absorber tube is 
enhanced at the presence of nanoparticles. The heat transfer 
coefficient is increased as the concentration of the nanoparticles 
in the base fluid is increased. For a given inlet Reynolds number, 
the enhancement in heat transfer coefficient caused by nanopar- 
ticles is decreased as the operational temperature of the absorber 
tube is increased. The maximum heat transfer coefficients were 
obtained at the left and right sides of the tube (@=310 and 
0=230).The results of this work for investigating the heat transfer 
variations of the assumed collector were validated with the other 
similar studies. 

Utilizing nanofluid in PTC has environmental benefit and 
causes improving heat transfer, consequently reducing the needed 
heat transfer area of the tubes and heat exchangers. Nevertheless, 
these advantages could be achieved with little change (in terms of 
materials, system design, and initial capital investment) on the 
entire PTC system. 


Appendix A 


Fluid properties: 


Por = —0.99854T + 1236 kg/m? (1) 
P; = 3850 kg/m? (2) 
Pef = 0 -HP + Ops (3) 


The correlation of effective density is taken of [33]: 
Hyp = 6.6710 ~7T* — 1.564107 °T? + 1.38%T? —5.541%10°T +8.848:104 uPa s 
(4) 


Hopp = (1 +2.50)flyp (5) 
The correlation of effective viscosity is taken of Einstein's [34]: 


Kop = —5.753496%107 !°T? — 1.875266%10~4T + 1.900210%*10~ ' W/m K 


(6) 

Ks = 5.5+34.5xe(— 0.0033«(T — 273) W/m K (7) 
2K, —2(1+ 8)? = 

Keff = Kp ee e (8) 


Ks +2kp+(1+ BY Ply —Ks) 


where / is the ratio of nanolayer thickness to the nanoparticles 
diameter which is considered to be 0.1 in Eq. (8) of this study [35]. 


Cpyp = 0.001708T + 1.107798 kJ/kg K (9) 
Cp, = 1.046 + 1.744107 T —2.79%104T? kJ/kg K (10) 
(CP) orp = A — P)CDpy + PCPs (11) 


The correlation of effective specific heat is taken of [33]. 
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Table A1 
Nanofluid (Al203) thermo-physical properties at 300, 400 and 500 K. 


Property T=300 K T=400 K 
0 1% 3% 5% 0 1% 
p (kg/m?) 936.5 965.6 1023.9 1082.1 836.6 
Cp (J/kg K) 1620 1612 1595 1578 1791 1782 
K (W/mK) 0.134 0.139 0.150 0.161 0.115 
u (mPas) 6.68 6.84 7.18 7.518 2.16 


In all of the correlations f, s and eff correspond to the based 
fluid, nanoparticle and nanofluid, respectively and @ is the volu- 
metric concentration of Al,O3 in the nano-fluid. The thermo- 
physical properties of the nanofluid were calculated via (3), (5), 
(8), (11). The governing equations for the fluid motion in a 
cylindrical coordinate system are described as follows: 

Continuity equation: 

10 


ror a2) 


lo 0 
TOCH 5p Petsr™Y) + Pep W) = 0 
Momentum equation: 

0 component: 


lo loa ð 1 
rage) + T greaf + a7 Pett WU) +7 Pefplv) 


1/fop\ 10a ou ð (Hepp oru) 
= (33) taalaa) a r or ) 


2Merf OV 
ji ai 99 PerrEBer Tw —T)sin 0 


a3) 


r component: 
12 12 0 1 2 
rage) + sea Hoz WV) ~~ plPnsl ) 
1/op 10 ov ð (nf rv) 
(2) tp erag) +3( r or 


Ung ou 
= 20 —P rf SP rp Tw —T) cos 0 


(14) 
Z component: 
12 10 ð 5 
TOLEN +7 r LAN + Pep W) 
op 10a ow\ 10 k ow 
az 72 00\ Hell ag) tr ar\ Melt ar 


Energy equation: 


(15) 


1a 192 0 
7 glen Hrga YD Hzer WT) 


10 ( Keff 5) 4 0 (: Keff z) 

~ 72 90 (CP) orp OO} ror \ (CP)eff OF 
where r, 0, and z are the coordinates, u, v, and w are the velocity 
components and p, T and g are the pressure, temperature and 
gravitational acceleration. p, k, y, and C, are the effective density, 
thermal conductivity, dynamic viscosity and specific heat capacity 


of the nano-fluid. 
The fully developed entrance length: 


1.22 < Lp =(4.4%Re!/®)xD < 1. 
6L: - 10xD = 0.7 m 


See Table A1. 
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